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ABSTRACT 


The object of the thesis is to study the role of flood plains 

on flood peak subsidence and time of occurrence of peak. A 

two-dimensional mathematical model for flood routing through a 

river with flood plains has been formulated. The model uses only 

continuity equation. By two-dimensional, it does not refer to the 

unsteady flow equations in two spatial dimensions, but rather to 

the physical situation in which river and flood plain cells form a 

two-dimensional network in the horizontal plane, but flow of water 

between two ceils is one- dimensional . The results of the present 

model are compared with the results obtained by liahapatra, K. , 

(1990) for one— dimensional model with both continuity equation and 

equation of motion. Laws of exchanges of flow between two ceils 

are governed either by river type link or weir type link depending 

on natural topography. River type link is based on Manning's 

equation, whereas weir type link follows the classical weir 

formulae. Since slope of channel is steep, bed slope has been* 

used rather than water surface slope in Manning's formuia- 

The parameters governing the problem have been identified 

through a dimensional analvsis. Values of B = (B - B )/B = 4, 8 

r s m m 

and = n^/n^ = 1,4 have been used and the other parameters have 

been kept constant. The computer program has ■'been developed in 

Fortran IV language and implemented on HP— 9000 and micro— vax. 

From the computer results for all the sets of data it is found 

that with = 1.0 model does not exhibit any flood peak 

subsidence as flood wave propagates down the river, whereas with 

n = 4.0 the model exhibits flood peak subsidence- After 

r 



comparing the results of the two dimensional model with the 

results of the one— dimensional model it is found that n has more 

r 

effect in two-dimensional model on flood peak subsidence whereas 
has more effect in one— dimensional model. The present study 
uses explicit method and some anomalies have occurred. It is 

recommended that the same problems should be studied using 


implicit method. 



ACKNOWLEDGEMENT 


-The author is highly thankful and sincerely grateful to his 
lesis supervisor Prof- S. Surya Rao for his valuable guidance, 
aspiration and encouragement during the course of this work. 

The author is extremely thankful to Dr. B.S. Murty of the 
apartment of Civil Engineering for giving him valuable 
jggestions and instructions during computational work. 

Sincere thanks are also extended to Dr. B. Dutta, Dr. T. 
angadhariah. Dr. V. Lakshminar ayana. Dr. B.S. Murty, Dr. S. 
amaseshan. Dr. K. Subramanya and Dr. S. Surya Rao of the 
epartment of Civil Engineering for the thought provoking courses 
3 has received from them during the M. Tech- Programme. 

The encouraging and unreserved help received from all his 
riends specially from Shri K.V. Jaykumar, Shri Ramesh Kumar 
rivastava, Shri Pranav Mahapatra, Shri D.K. Purvar, Shri 
ivekanand Singh, Shri Y-J-J. Frasad, Shri R-R- Rakesh, Shri G-P. 
ingh and Shri K-R-C. Rao Kalaga during the various stages of the 
hesis preparation is thankfully acknowledged - 

Last but not the least the author wishes to place on record 
is sincere thanks for the kind assistance and cooperation 
rovided by Shri S.K. Jlishra, Shri Suman Kumar and Shri Shambhu 
ha during the entire degree programme. 


Sunil Kumar Choubey 


UST OF RGURES 


Fiour 

e_ 



Pane 

2. 1 


Continuity equation for a cell- 


37 

2.2 


Weir type link. 


38 

2.3 


Arrangement of cell groups. 


39 

3. 1 


Channel cross-section with flood piains. 


40 

rr '-y 

^ m ^ 


Inf low hydrograph .. 


40 

- 3.3 


Topological breakdown of computat ional cells 

- 

41 

3.4 


Illustration of oscillations in weir flow between 




two cells. 


42 

^ • 5 

(a) 

Stage hydrographs in main river for B^=4.0, 

n =1-0. 
r 

43 


Cb) 

Stage hydrographs in flood plain for 

B = 4 . 0 , n = 1 . 0 . 
r ’ r 


44 


(c ) 

Stage hydrographs for one— dimensional model 

B =4.0, n= 1-0- 
r r 

for 

45 

M 6i 

(a) 

Stage hydrographs in main river for B^=4.0, 

n =4 - 0 - 
r 

46 


(b ) 

Stage hydrographs in flood plain for 

B = 4.0, n = 4.0. 
r r 


47 


(c) 

Stage hydrographs for one— dimensional model 

for 




B = 4.0, n = 4-0. 
r r 


48 

3.7 

(a) 

Stage hydrographs in main river for B^=8.0, 

n =1-0- 
r 

49 


(b) 

Stage hydrographs in flood plain for 

B =8-0, n = 1-0- 
r ’ r 


50 


<c> 

Stage hydrographs for one— dimensional model 

B =8-0, n=l -O- 
r ’ r 

for 

51 

3. S 

(a) 

Stage hydrographs in main river for B^=8-0, 

n =4.0. 
r 

52 


(b ) 

Stage hydrographs in flood plain for 

B'=8-0, n =4-0- 
r ’ r 


53 


(c) 

Stage hydrographs for one— dimensional model 

B = 8.0, n = 4-0- 
r ’ P 

for 

54 

3-9 

(a) 

Normalized time of peak depth vs- normalized 




distance along the main river. 


55 


(b) 

Normalized time of peak depth vs. normalized 



distance along the flood plain. 


56 



3.10 (a) 

Normalized peak depth vs normalized distance 
along the main channel. 

57 

(b) 

Normalized peak depth vs normalized distance 
along the flood plain. 

58 

<c) 

Normalized peak depth vs. normalized distance 
for one— dimensional model. 

59 



(vin) 


UST OF SVMBOLS 

A. . = Area of the flow cross section between cells i and k. 

1 ,k 

AS^ = Water surface area in horizontal plane of the i—th 

' cell - 

AAS^ = Increment of water surface area due to change in 

water level- 

B = (B -B )/B . 

r s m m 

B = Width of main channel. 

m 

B^ = Width of channel including flood plains. 

Cdw = Coefficient of weir- 

Cdo = Coefficient of orifice. 

dt = Time increment. 

dx = Distance increment- 

dZ. = Water level increment. 

1 

g = Acceleration due to gravity- 

i,k = Subscripts indicating cells- 

K = Conveyance of the cross section - 

= Manning's n of main channel . 

= ManninQ's n of flood plains. 

"r ^ n^/n^. 

P. = Rainfall intensity. 

Q. . = Flow exchange between cells i and k. 

1 ,K 

R. , = The hydraulic radius of the flow cross section between 

cells i and k. 

S = Water surface slope - 

S 
o 


= Bed slope. 



t 


t 

p 

V. 


eak 


AV. 

1 

X 

y 

y' 


z. 

□ 

z, 

X 


7 


W 

2w 


o 

Zwo 


P 


Time under consideration . 

The time to peak stage of inflow stage hydrograph- 
Time base of in flow hydrograph. 

Volume of water stared in cell i. 

Increment to volume of water stored in ceil i. 
Distance along the channel. 

Depth below bottom of flood plain. 

Steady state depth. 

Water level in the river cells above river bed level 
section under considerat ion . 

Water level in the flood plain cells above river 
level at that section. 

Bed level from datum. 

Water level in cell i. 

Elsvation of water surface. 

Steady state stage. 

peak stage of inflow hydrograph. 


at 

bed 



TABLE OF CONTENTS 


CERTIFICATE 
ABSTRACT 
ACKNOWLEE>GEMENT 
LIST OF FIGURES 

LIST OF SYMBOLS 

CHAPTER I INTRODUCTION 

1 . 1 General 

1.2 Review of Literature 

1.3 Scope ot Present Study 

1.4 Closure 


Page 

(i) 

(ii) 

(iV) 
Ctxii ) 


1 

4 

6 

7 


CHAPTER II GOVERNING EQUATIONS AND METHODS OF SOLUTION 


2-1 Introduction 8 

2.2 Assumptions 8 

2-3 Darivation of Equations 9 

2.3.1 Continuity Equation for a Ceil 9 

2-3-2 Laws of Discharge Between Ceils 11 

2.4 ^4umerical Solution of Governing Equations 14 

2.4.1 Numerical Formulation 15 

2.4.2 Numerical Solution of Difference Equations 17 

2.5 Accuracy and Stability 19 


CHAPTER IX-I TWO-DIMENSIONAL MATHEMATICAL MODEL FOR FLOOD 
ROUTING THROUGH A RIVER WITH FLOOD PLAINS 


3-1 General 

3.2 Channel Geometry 

3.3 Initial Condition 

3.4 Boundary Condition 
3.4-1 Up Stream 

3-4.2 Down Stream 

3.5 Model Building 

3-6 New Depth Computation 
3-7 Stability Condition 
3.8 Computer Program Structure 
3-9 Computational Difficulties 
3.9.1 Weir Oscillations 

3.10 Selection of Input Data 

3.11 Computations 

3-12 Results and Discussions 


20 

21 

21 


24 

24 

25 
28 
28 
29 
31 



CHAPTER IV CONCLUSIONS AND RECOMMENDATIONS 


4.1 Conclusions 35 

4.2 Recommendations 36 

FIGURES 37 

REFERENCES 60 

APPENDIX A 61 



1 


CHAPTER I 

INTRODUCTION 


1 . 1 GENERAL 

Mathematical modelling provides a tool by means of which one 

can : 

(i) Study and gain an understanding of hydraulic flow 
phenomena- 

(ii) Select and design sound engineering projects. 

(iii) Predict extreme situations so as to be able to provide 
advance warning of their occurrence and importance, and 

Civ) For proper management of flood plain area. 

The prediction of exceptioTial natural events and their 
consequences implies the use of mathematical modeliing because 
past observations of such events are so often non-existent. The 
capability of a properly constructed mathematical model to perform 
well depends on the fact that it uses physically correct, 
determin ist ic - hydraul ic equations whose validity is not limited to 
the known flow events used to adjust their empirical coefficients. 

Flood Routing is basically determination of stage hydrograph 
at any section in the river if the hydrograph (stage or discharge) 
is given at an upstream location. There are two most important 
parameters in analysing the flood (i) Peak depth and (ii) time of 



occurrence of peak depth. There are several methods available far 
predicting these parameters for any given inflow hydrograph and 
channel parameters. These models may be broadly classified under 
two headings: 

(i) One-Dimensional Model and 

(ii) Two-Dimensional Model 

Truly straight channels in which the flow can be considered 
to be strictly one— d imensional are seldom observed in nature. It 
is more often the case that over bank flows go their own way, 
filling up the flood plain in a manner dictated by local 
topography, sometimes never returning to the main channel during 
the falling flood. 

By two dimensional we do not refer to the unsteady flow 
equation in two spatial dimensions, but rather to the physical 
situation in' which channel and storage cells form a 
two-dimensional network in the horizontal plane, but flow of water 
between two cells is one dimensional. 

Mathematical modelling of unsteady flow in river is much more 
complicated than in artificial channels because of possible 
inundation of flood plain and consequently flow is not 
one-dimensional any more. The construction of the mathematical 
model of a flooded area sets up a series of p rob lems concern ing 
both the representation of the physical reality and the numerical 
methods. Many times it is not passible to predict the extent of 
inundation itself. Even with well defined flood plains the flow 
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pattern is not simple and interaction between the main river and 
the flood plain cells affect the flood wave propagation. 

During the rising flood, the positive wave advancing along 
the river will have higher celerities and higher stages along its 
length. Thus there will be a transverse slope which will fill up 
the inundated plain with water. During the falling flood, the 
reverse situation will be observed: the main channel will try to 
drain first and the transverse slope will then cause water to flow 
from the inundated flood plain back to the main channel. But one 
dimensional modelling as described by the St. Venant equation 
written for nan-prismatic channels implies that the transverse 
slope is nil. The lateral flow influences the celerity and 
consequently the time of occurrence of peaks (Cunge, 1975). The 
flood is also attenuated because of this. Therefore any accurate 
model should take into account these two-dimensional effects. If 
in the most situations purely one— dimensional models are used, it 
is because two-dimensional ones have historically been either too 
difficult or too expensive to construct and apply. It is fortunate 
that most often the overall results of these models do not 
seriously differ from reality, although they cannot always 
reproduce complicated details of flow which are ail but 
one- dimensional . The present work is an endeavour to develop a 
two-dimensional mathematical model to study the role of flood 
plains on the flood wave propagation in a river with flood plains 
and also to compare it with one dimensional model. 
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1.2 REVIEW OF LITERATURE 

Mathematical modelling of unsteady flow for ideal condition 
which uses the complete equations is very common and widely 
covered in several text books (Gunge, et.al., 1980, Chaudhary, 
1987). However, technique finds its origin in the 19th Century 
work of de St. Venant and Boussinesq who formulated the unsteady 
flow equations, and in the work of Massau who in 1889 published 
some early attempts to solve those equations. Most of the methods 
are applicable to only ideal situations like simple channel. V/ery 
few models have been developed to meet practical situations. 
Important theoretical concepts were established in first half of 
this century, but the first engineering applications of this 
principle to the natural river conditions awaited the development 
of digital computers. Isaacson, E. , Stoker, J., and Troesch, 
B.A., (1954) constructed and ran a mathematical model of portion 
of the Ohio and Mississippi rivers by means of finite difference 
method on the digital computer, UNIVAC. Special attention was 
given to problems resulting from junction of the Ohio and 
Mississippi rivers. Following that pioneering effort the use of 
mathematical modelling in natural river began to accelerate. 

Gunge, J.A., (1966) presented a paper discussing the dif f ic^li ti es 
encountered in the development of mathematical models for the 
reach of Oraison canal. He compared the mathematical model with 
scale model and comparative results were presented. He concluded 
that mathematical model results are quite comparable with those 
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given by a distorted scale model and use of mathematical model in 
preference to distorted scale models is justified. 

Later on Cunge, J.A., Lorgere, H. , and Preissmann, A., (1966) 

constructed the UNESCO mathematical model for the control of 


floods of 

the 

river Mekong. 

Zanobetti, D 

. , Lorgere , 

H., 

Preissmann , 

A. , 

and Cunge, 

J .A. 

, (1968) were 

involved in 

the 

deve lopment 

and 

app 1 icat ion 

of SOGREAH'S CAR IMA 

system for 

the 

simulation 

of unsteady flow 

in 

Mekong river 

and Mekong 

flood 


plains. In 1970, Zanobetti , D., Lorgere, Henri, Preissmann, 

Alexandre and Cunge, J.A., constructed and used the model for 
lower basin of the Mekong river. The model consisted of a system 
of ordinary first order differential equations each representing 
water level as a function of time in one of the some 300 meshes in 
to which whole area was divided. Exchange relations between ceils 
were based either on Saint Venant's dynamic equations or on the 
Weir type exchange laws, inertia terms being ignored in both cases 
due to the very small velocity of flood propagation. Cunge, J.A., 
(1971) presented a paper in proceedings, 14th Congress, lAHR, 
Paris discussing the representation of flood plains in the 
mathematical models in rivers. Grijsen and Vreugdenhil, (1976) 
develop-ed a mathematical model of the Rharb plain in liorocco. 
(Basically mathematical formulation is as in one dimensional 
case). In two-dimensional model they considered inter connection 
of one cell with its neighbouring cells. In the case of the Oused 
Sebou and the Rharb plain, the network approach was used and full 
dynamic equations including inertial effects were used. But they 
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realised that if there is gradual variation of discharge or stage, 
then inertia does not play an ^ important part. Cunge also 
concluded similar idea about neglecting inertial effect. 

Tingsanchali and Ackerman (1976) derived equations for 
compound channel flow which include both dynamic and storage 
effects of the over bank flow. They applied these equations to 
study the floods in a natural river and concluded that higher peak 
stages and lower peak discharges are obtained when only storage 
effects of flood plains is considered. Later Tingsanchali and Lai 
(1988) using the same equations obtained semiemperical exponential 
equations to describe the subsidence of flood peak discharge in 
channel during over bank flow periods. 

Finally, Hahapatra, K. , (1990) routed the flood through 

compound channel using inertial effect and concluded that one can 
neglect inertial effect for natural river and flood plains. He 
reccmmerHied for development of two-dimensional mathematical model 
of compound channel ignoring the inertial effects. 

1.3 SCOPE OF PRESENT STUDY 

1. To present a two-dimensional mathematical model of a river 
with flood plains. Inertial effect has been ignored from 
dynamic equation. An explicit finite difference method is 
used to solve governing equation numerically. 

2. To study the role of channel parameters on the flood wave 
propagation. 
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3. To do comparative study with one-dimensional model with 

inertial effect developed by Mahapatra, K., (1990) » 

1 . 4 CLOSURE 

The following chapters present the development of 
two-dimensional mathematical model of river with flood plains- 
Chapter II presents governing equations and method of solution- 
Assumptions made are also given in Chapter II- Chapter III 
explains the application of model to a river with flood plains 
results and discussion- Chapter IV includes conclusion and 
recommendations for future work- The computer programme has been 
included in Appendix A- 



CHAPTER II 


GOVERNING EQUATIONS AND METHOD OF SOLUTION 

2.1 I NTRODUCTI ON 

In this chapter the governing equations for unsteady flaw 
without inertial effect have been presented. In addition to this 
method of numerical solution to these equations has been given. 

The arrangement in this chapter is as fallows. In Section 

2.2 assumptions made in the derivation of the aquations are 

discussed. In Section 2.3 continuity equation for a ceil and laws 
of discharge between the cells are presented- In Section 2.4 

numerical solution of governing equations are presented and in 
Section 2-5 accuracy and stability are discussed - 

a^a ASSUMPTIONS 

The assumptions made for deriving the governing equations are 
as follows: 

(i) There are gradual variations in flow characterist ics , so 
that one can neglect inertial effect. 
ill) The volume of water stored in cell i is directly related 
to the level in this cell i.e., 

V. = ¥ <Z.) 

11 
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(iii) The discharge Q. between two adjacent cells i and k at the 

1 , k 

Given time t = n A t is a function of the levels and 
^ n i k 

only, i • e . 

= 0 ( Z-, Z- ) 

(iv) The discharge between two calls Q, .. does not vary during 

1 , K 

the period of time At. 

(v) The slope, 0 of the channel bed slope is small so that sin & 

tan ^ 3 and cos S ^ 1, 

(vi) The resistance coefficient for unsteady flow is same as for 
steady and uniform flow. 

(vii) The velocity distribution in a vertical plane is uniform. 

2.3 DERIVATION OF EQUATIONS 

In this section continuity equation for a cell and law of 
discharge between cells are derived. The derivation closely 
follows the derivation given by Cunge J.A. (Chapter 17, K. Mahmood 
and V. Yevjevich). 

2* 3* 1 Coritinuity Equation for a Ceil 

Let us consider a ceil as shown in Fig. 2.1. 

Suppose that to every cell i corresponds a characterist ic 

water level assumed at the center of the ceil. It is also 

assumed that the water surface is horizontal within the cell 

boundaries and its elevation above the datum is Z . . 

1 

A constant time step At is assumed. Then, for any given time 
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A constant time step At is assumed. Then, for any given time 
t = nAt the water in cell i is Z..(tn) and the corresponding water 

n _ 1 I- =F 

surface area (the surface ABCD) is equal to AS. (t ) the water 

i n 

surface area in the horizontal plane of the cell. 


At the next time step t . = t + At the water level is Z. 

n+1 n 1 

<t ) and water surface area is As.(t .) = AS.(t ) + AAS . - The 
n+l 1 n-*-l in i 

increment of water level may be due to runoff P^{t) resulting from 
rainfall P^(t) on the cell surface during the time At and to 
discharge Q. . (Z. ,Z. ) from the cells K adjacent to the cell i. 

X f K X K 

The increase in water volume stored in the cell i during the time 

At from discharge conditions is 

t -^-At t -J-At 

n n 


AV. 

1 


J P <t) dt + E J Q 
t t 


n n 

where, the suffix k represents numbers 
cell i. From geometric considerat ion , 
written as below 


i,k '^i» '2.1) 

of all cells adjacent to 
the above relation can be 


Z. (t +At) 

1 n 

AV. = f AS. (Z.) dZ. (2.2) 

1 J 111 


Z . (t ) 
1 n 


Suppose that the surface area AS. does not change between the 

ZiAS. ^ 


levels Z. and Z. + AZ . i.e, 
1 11 


( 


AS 


- « 1 ) 


From Fig. 2.1, 


AV. = V 


ABCDA'B'C'D* 


= V 


ABCDA"B"C"D" 


= AS. (Z. ) AZ 
11 ] 
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This relationship replaces integral sign in equation <2.2) 
and dZ^ is replaced by 

Equating equations (2-1) and (2.2) 


AS. (Z . ) AZ . = P . (T) . At + At 2 Q . , £Z . (r) , Z, (T)> 

11 1 1 j^i,ki’K 

. where , 

n.At < T ± (n+1) At 


If AZ -* 0, and At ^ O. 

Equation (2.3) can be written under the differential form 
dZ . 

J Qi.k 

k 


This is the continuity equation for cell i- 


(2-3) 


(2.4) 


Similar equation can be written for each cell. Thus for N 
cells there would be N ordinary differential equations for N 
unknown functions 2^ of the independent variable t. 

2.3.2 Laws of Discharge between the Cells 

Generally two types of connections between ceils are assumed - 
These are: 

Ci) River type link 
(ii) Weir type link 

River type links represent the exchanges between two cells 
when there is no local obstacle to the flow and a mean resistance 
co—efficient for a given cross section of the flow can be used- 


! 
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Weir type links represent exchanges between cells where roads 
or dikes form -a local obstacle that can be represented by singular 
head loss between two cells. 


2. 3. 2. 1 River type link: 

In the case of river type link Manning's formula is used. 


Q. 


i ,k 


n 


A. , 
1 ,k 


/3 
1 , k 


.1/: 


(2.5) 


where , 

n = Manning's roughness coefficient 

A. . = the area of the flow cross section between cells i 

i,k 

and k 

R. = the hydraulic radius of the flow cross section 

1 , k 

between cells i and k, and 
S = water surface slope- 


In equation <2-5 ) parameters n, A. and R. are function of 

X f rC X ^ K 

water level 2' , in the flow section between cell i and k. We can 

i,k 

write 


- . A. , . = K(Z. , : 

n i,k i,k x,k 


where , 


Z . . = d Z . + ( 1-a) Z, 

1 , k 1 k 

Z. . is weighted level between levels in two cells i and k and 

1 , K 

KCZ. ,.) is called as conveyance. 

i , K. 


a is called weighting coefficient and it is constant for a 
given pair of cells. For equal sizes of cells, for given pair 
is equal to 0.5. For bigger size of cells weighting coefficient 
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more than 0.5 is provided for that cell. 

The water surface slope S in equation (2.5) is given as 



where , 

Ax is the fixed distance between the centers of cell i and k. 
Let us define a function 




kCZ. , ) 

-TST 


n*TSx 


where, 

k is conveyance of the cross section. 

The discharge formula can be written as 

ai,k = - Zii . -1(2, - Z-) 

Sign will depend upon the convention adopted with regard to 
the direction of flow- The function must be established 

fC 1 

first based an the channel parameters. 

2. 3. 2. 2 Weir type link 

In this case the classical discharge formulae for weirs are 

used . 

There are two, types of flow, namely submerged flow , and free 
flow depending on following criterion: 


Considering Fig. 2.2<a) 

if w < ^ (w+h), the flow is free flow and discharge is given by 


Q 


i,k 



3/2 


42g . Cdw . dx . h 


( 2 . 6 ) 
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if w > 2 tw+h), the flow is submerged flow and discharge is given 

by - _ 

Q. = J.42Q.Cdw . dx . + 4^-CdcT . dx . w (2.7) 

X f K! 1.^ 

where , 

dx = length of Weir (each cell) 
h,w = as shows in Fig. 2.2(a) 

Cdw and Cdo are coefficient of Weir and coefficient of orifice 
and given by Rouse as follows 

Cdw = 0.611 + 0.075 h/w 
Cdo = 0.611 - 0.175 (w/h+w) 

a. 4 NUMERICAL SOLUTION OF eOVERNING EQUATIONS 
Introduction 

With the advent of digital computers, more emphasis is being 
given on numerical solution methods. However, there have been 
some analog solutions. While using the numerical scheme, one 
should bear in mind the convergence, stability, efficiency and 
accuracy of the scheme before selecting it for a particular 
application. The governing equations have become very simple 
after dropping momentum equations. It has been left only with 
continuity equations and they- are a set of ordinary differential 
equations. So, there is no need to use any special numerical 
technique. The set of differential equations can be solved by 
finite difference method either by explicit method or implicit 
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method dspendinQ on the assumption made in variation of Q- . 

1 , K 

during the period of time At . ’ 

2. 4. 1 Numerical formulation. 

There are two ways of writing last term of equation (2.3) 
based on different assumption e.g. (i) Explicit method <ii) 

Implicit method. 


2. 4. 1 - 1 Explicit method 

It may be assumed that the discharge Q. does not vary 

1 , K 

during the period of the time At, i.e. 

T = n . At 

and 


Q. IZ.Ct), Z, <t)> = Q.'t^in.At), Z, (n.At)> 

Xi|rCX rC X^K K 

In this case equation (2.3) becomes 


Q 


n 

i,k 


AZ. . AS. (Z") = P'l' . At + At 2 , (2.3) 

1 11 1 If 

It means, if the levels in the neighboaring ceils are known 
at time tn = n. At- It gives an explicit solution for AZ^, because 
all terms on right hand side of equation are known, then the 
resulting equation can be solved for each ceil at a time. 


Consequently, the water level at next time step (n+l)At will 


be 


z""^^ = z" + AZ. 

1 1 X 
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2. 4. 1.2 Implicit method 

It may be assumed that the discharge Q. , (Z . (t) ,Z. (t) ) is an 

X ^ K X rv 

intermediate discharge between Q? . and . In other words, 

X a K X a tC 


nAt < T < (n+l ) At 



Since, the inter cell exchange laws (see equation river link, 
Weir link) are nonlinear and their substitutioh into equation 
(2.10) mould further complicate the problem. For each cell there 
would be one nonlinear algebraic equation like (2-10). Since the 
model will have several hundred cells, consequently a large number 
of nonlinear equations will be obtained- It would not be a easy 
job to solve a large number of algebraic equations- In order to 
make the method of solution easier, the difference equations are 
linearized using the Taylor's series expansion and neglecting 
higher order terms with an assumption that the water level 
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Thus, a system of linear equations for the is obtained. 

1 

After substituting the value of in to" equation (2.10) 




dl. 


^i,k 

&L, 


( 2 . 12 ) 


Since the above equation contains two unknowns AZ^ and 
one cannot get AZ^ directly from it. Similar equation can be 
written for each cell and system of linear algebraic equations is 
obtained. The system of finite difference used in Eq . (2.12) is 

the implicit type depending on the value of weighing coefficient & 
to express the derivatives. This is very important as far as the 
numerical solution is concerned. 


By putting 8=1 equation (2.12) becomes 


AZ - 

as” = p. 
1 At 1 


^ ^^i k ^ 

i ? K 1 . 


AZ. + Z 
^ k 


i ,k 


(2.13) 


To find the values of AZ . and AZ. the algebraic equations for 
the entire system of ceils have to be solved simultaneously- 
Impiicit finite difference schemes usually do not have to satisfy 
time constraints and therefore, longer computation time steps can 
be used. They involve matrix inversion and require large number 
of calculations, somewhat masking the advantage gained by using 
longer time steps- 


2. 2 Miimerlcal Solution, of* Difference Equations 

seen in the Section 2-4-1, the continuity equation (2-3) 
can be written in two ways, either in explicit way or in implicit 
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way depending on assufnption made for Q. For solving these two 

1 , K 

equations there are cprrespondingly two methods namely (i) 
Explicit method, (ii) Implicit method- 


2. 4. 2. 1 Explicit method 

Rewriting the explicit type equation <2. 8) 

AZ. AS. (Z'^) = At + At Z q" . 

Ill 1 j^i,k 


AS = P" + Z q" . 
dt 1 1 , i,k 

k 


As it is clear from the above equation, AZ^ can be obtained 
without any need of other equations. Its solution is nothing but 
a solution of ordinary differential equation. So any numerical 
method for solving O.D.E. is sufficient. Thus explicit method 
simplifies the problem drastically but it has got a great demerit 
relating to stability. The price paid for this simplification is 
the limited choice of time step At. So, a considerable amount of 
computer time would be necessary. Thus for long period analysis 
implicit method is preferred in order to obtain a free choice of 
time step At. 
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The cells 
in Fig. 2.3 in 
group directly 




.k 


are arranged in 
such a way that 
above and below 


certain number of groups as shown 
\ 

each group is related only to the 
it. 


Equation similar to (2-14) can be writt:en for t^le central 
group iy the preceding group j — 1 and the folloufing group i + 1. For 
details one may go through Chapter 17 of "Unsteady Flow in 
Channels", Vol . II, edited by K. flahmood and V, Yevjevich- 
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CHAPTER III 

TWO-DIMENSIONAL MATHEMATICAL MODEL FOR FLOOD ROUTING THROUGH 

A RIVER WITH FLOOD PLAINS 


3. 1 GENERAL 

Unsteady flaw in a river with flood plains is usually two- 
dimensional if the flow depth in the flood plain is shallow. By 
two dimensionality, we refer to the physical situation in which 
cells form a two-dimensional network and the exchanges of water 
between cells are purely one dimensional. For exact analysis, 
two dimensional modelling is required. Due to very limited time 
a channel with well defined flood plain has been chosen for 
study. 

The aim of present work is to study the effect of variation 
of parameters of flood plain character izing the geometry and flow 
condition during over bank flow periods- In addition to this, 
the results of this model are compared with the one dimensional 
model developed by Hahapatra, K. (1990) in which he has taken 
inertia terms into considerat ion . Inertia terms have been 
neglected assuming that variation in flow characteristics is 
gradual which is more likely to be practical situation. 
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3.2 CHANNEL GEOMETRY 

As shown in Fig. 3.1-, a symmetric prismatic channel with 
flood plains has been chosen for the study, primarily with a view 
to compare the results of the two-dimensional model with the 
results of the one— d imensional model studied by liahapatra. The 
depth below the flood plain bottom is 5 m. There is no 

transverse bed slope in the main channel as well as on flood 
plain. The longitudinal bed slope for main channel and flood 
plain are same. The total length of channel reach is 10 km but 
only first 8 km has been studied to prevent any back water effect 
or down stream boundary effect. 

In Fig . 3.1, 

= width of main channel 

= width of channel including flood plains 
= depth below flood plains 
= depth of flow 
= elevation of water surface 
= bed level from datum. 

3.3 INITIAL CONDITION 

The depth of flow at the center of each cell is defined at a 
given time (t = 0). Initially steady uniform flow is assumed 
with water level slightly more than the level of bottom of flood 
plains. The exchanges of flow have been computed based on river 

j 

type link or weir type link depending on physical topography. 
River type link is governed by Manning's equation. The weir type 


B 

0) 

B 

s 

i 

y 

w 

"“b 
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link is governed by Weir formula, given by Rouse. The channel 
bottom slope is assumed to be constant along whole reach and 
throughout the duration of interest. 


3.4 BOUNDARY CONDITION 


3. 4r. 1 up stream 

For the present study a symmetric triangular shape of stage 
hydrograph of 1 hour has been considered at the upstream boundary 
cells- The shape of hydrograph is shown in Fig. 3.2, 
where , 


wo 


wop 


peak 


t 

t. 


steady state stage 
peak stage 

the time to peak stage 
time under consideration 
time base of hydrograph 


The water level in all cells at upstream boundary is same at any 
time. It means there is no transverse water surface slope in 
upstream cells. The new water surface level is calculated 
directly from stage hydrograph for new time. 


3- 4- 2 Down Stream 

The new water level at down stream cell is computed by 
assuming that the change in water level in down stream ceil is 
same as that of just preceding cell. 

3-5 MODEL BUILDING 

The division of an area like river with flood plains into 
individual cells, same of which are portions of main river and 
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other of the flood plain must be done carefully. The breakdown 
of given area into cells is based on natural topography. Since 

in present case there are well defined flood plains and main 
river, division of an area into cells is not a difficult task. 
The cells have been divided into four categories. They are main 
river cells, ne;<t to the river flood plain cells. Intermediate 
flood plain cells, and Side boundary cells. After dividing the 
whole area into cells- the picture looks like Fig. 3.3. 

As shown in Fig. 3.3 each cell is linked with neighbouring 
cells with a link either river type or weir type depending on 
geographical topography. Since there is abrupt change in bed 
level between the river cells and just next to river cell, so 
they are connected with Weir type link because there would be 
singular head loss. Remaining links are river type links. River 
cells carry high flow compared to flood plains cells. The cell 
centers are their geometrical centers. However it is not 
necessary condition and it should be defined in such a way that 
the directions of flow are correctly allowed for any time during 
flood. The sizes of the flood plain cells and river ceils are 
same and rectangle in shape of dimensions 500 m x 50 m. The bed 
level and water level have been assigned at the center of cell 
and assumed horizontal throughout the cell. The water surface 
area of ceil is not a function of water level but it is constant. 
The manning's n is not a function of stage but constant. The 
side wall has same value of manning's n as its neighbouring bed 
surface. In river type link equation bed slope has been used for 
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flaw in the longitudinal direction whereas water surface slope 
has been used for flow in the transverse direction. 

3.6 NEW DEPTH COMPUTATION 

After neglecting the rainfall terms in equation (2.8>, it 
becomes 


AZ . 
1 

At 


AS. 


= I 


(3.1 ) 


Above equation is an algebraic equation and the value of AZ^ 
is found explicitly from given relation. 


For any interior cell, all the discharges connected to that 
cell with neighbouring cells are found based on water level of 
calls at known time step and the type of link. These discharges 
are algebraically added which gives the value of term on right 


hand side of equation 

(3.1). Since AS . 

1 

and 

At are 

known , 

the 

value 

of is found 

very easily with 

the 

following 

relation 



AZ. = Z q” , 
1 k ^ 

At 

AS . 

1 






The same procedure is 

repeated with all 

the 

cells. 

If 

all 

the 

cells 

are exhausted. 

the new values 

are 

reset 

and 

time 

is 


increased by At. Same procedure is repeated for all the cells 
till total time of computation. 

3.7 STABILITY CONDITION 

Since explicit method has been used for solution of the 
difference equation, it is very much Susceptible to unstability. 
The finite difference scheme is said to be stable if small 



25 


numerical error introduced at the beginning are not magnified 
during successive application of the difference equations and the 
errors at beginning do not grow so large that valid- part of the 
solution is obscured. So in explicit method of solution one will 
have very limited choice of selection of At and at the same time 
one should remain very much cautious while selecting the time 
interval At. 

3.8 COMPUTER PROGRAM STRUCTURE 

Explicit method for solution has been used for the present 
study. Runge Kutta method of 4th order has been used for the 
solving the ordinary differential equation. The flow chart is 
given later and the general program is given in Appendix A. 

The main program calculates the initial conditions, boundary 
conditions and solves the ordinary differential equation by Runge 
Kutta' s method. As mentioned earlier, there are four kinds of 
cells, namely, river cell, next to river flood plains cell, 
intermediate flood plains cell and side cell. Separate 
subroutines have been formulated for the calculation of exchanges 
of flow between cells for each kind of cell. The constant time 
step At of 2 sec. has been used. However, results are printed at 
intervals of 120 sec. The stage hydrographs are obtained at 2 km- 
intervals in main river as well as on flood plains. The peak 
stage and time of occurrence for those cells also have been 
obtained. Flow chart shows the algorithm used in the program and 
gives clear picture about it. 



Read and write 
channel parameters 


Set initial conditions 


Set the u/S boundary 
condition to find the 
water level in u/S celli 


Unsteady calculation at 
interior nodes starts 


Identify the cell type 


Call corresponding subroutine 
and find and 

in applying R.K. method and 
finally find dzw 


Calculate new water level 

and water depth at the 

end of that time level as 

Zw = Zw + dzw 
n 

Y = y + dZw 
n ^ 
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3.9 COMPUTATIONAL DIFFICULTIES 

In this section, computational difficulties, which were 
faced during developing the model and its remedies are discussed. 
The concern here is not with the formal mathematical behaviour of 
various numerical schemes but with the numerical treatment of 
particular situations related to the fact that we try to solve 
continuous flow system with a discrete model. One should bear in 
mind the fact the difficulties such as described below may not 
appear at all in certain modelling systems and on the contrary, 
other troubles may occur in such systems which are not mentioned 
here . 

3.9.1 Weir Oscillation 

In certain situations and for certain modelling systems 
which do not use an iterative solution procedure, the direction 
of flow over the Weir at a given time step can be in opposite 
sense from what it should be, and the persistence of the anomaly 
can destroy the calculation. 

Fig. 3.4(a), (b) and (c) represent three consecutive time 
steps. In Fig. 3.4<a) flood overflows arrive first in cell A, so 

that its level rises above at the end of the time step. During 
the ne;<t time step, the program which had not used an iterative 
procedure calculates a discharge from A to B. But due to 
tremendous amount of water exchange or small surface area in B, 
its level jumps to the position shown in Fig. 3.4(b) at the end 
of the time step. Of course , iterative procedures would correct 
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it bat if they are not used then at the beginnirog of third time 
step (b), the situation is anomalous in that the flow occurs in 
reverse direction i-e- from ceil B to cell A and at the end of 
time step <C) again the water level in cell A goes higher than 
water level in cell B. This type of oscillation generally 
diminishes especially if each of two cells communicates with 
others- But if the cells are isolated from other, the 
oscillation can become unstable and destroy the calculation- 

One solution to this kind of problem is to use iterative 
procedures to avoid oscillation, but it may take considerable 
computer time- Usually it is less costly to accept some 
inconsistency like the above oscillations provided that volume 
continuity is net violated beyond reasonable limits, 

3-9,2 Dif fer ment Alter- native Methods Attempted and ProbleiT^s 
Associated with Them 

In the beginning the problem has attempted taking water 
surface slope in the Hanning's equation for river links in the 
longitudinal direction as suggested by Cunge but the solution was 
found to be highly unstable- However^ Cunge did not mention in 
his work about such kind of difficulties. Hence, instead of 
water surface slope, bed slope was used in the Manning's equation 
for river links in the longitudinal direction, , 
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A number of different alternatives far upstream boundary 
condition were tried includinQ the discharge inflow hydrograph 
instead of stage inflow hydrograph. In the case of discharge 
hydrograph, compound channel approach was used to determine the 
depths in the river and the flood banks at each time step. 
Unnecessarily lot of iterations were required for binding out the 
depths in the upstream cells, consequently consuming more 
computer time without improving the results. To make the 
computat ion al procedure simpler, stage hydrograph was taken as 
upstream boundary condition. 

3.10 SELECTION OF INPUT-DATA 

A symmetric prismatic channel with flood plains as shown in 
Fig. 3.1 has been considered for the study- The length of the 
channel has been taken 10 km but only first S km has been studied 
to avoid the effect of down stream boundary. According to Fig. 
3.1, the main channel width B is equal to 50 m and the bank full 
depth is equal to 5 m. the manning's n of main channel has been 
taken as 0.03. The manning's n of the flood plain is varied from 
0.03 to 0.12. The slopes of the main channel and the flood 
plains are same and equal to 0.0005. the initial flew depth in 
main channel is 5-5 m and In flood plains 0.5 m. The datum line 
is assumed to be horizontal and it is lO m below the upstream 
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main river cell* The duration of stage hydrograph is 1 hr. and 
total time of computation is 8 hr.-- 

The relationship between the dependent variables and the 
independent variables of the problem under study are as follows: 




= f 


■f "T" ^ ^ 9 t . B . 

1,2,^ m 


B 


n 


1 ’ 


n^ ) 


^ I > 

peak 


where , 


y 


y 


t 

B 

m 

B^ 

n. 


'peak 


B , B 
m s 


the flood 
hydrograph . 


water level in the river cells above river bed level 
at section under consideration. 

water level in the flood plain cells above river bed 
level at that section. 

peak depth over a cell under consideration. 

distance of call from u/s boundary. 

time of computation 

width of main channel. 

total width of channel. 

manning's n of main river- 

manning's n of flood plains. 

- time up to peak stage of inflow hydrograph. 

, n. and n.^ are characteristics for the channel and 

banks, t , is characteristics of the inflow 

peak 


Dimensional analysis of the above listed variables will lead 
to the following functional relationships between the 


dimensionless groups: 
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-2- = ^ ^ B , n ) 

^0 ^0 ^peak *' ^ 


-2 = f ( 2i_ B , n ) 

8 y^’ r’ r 


peak 


where , 
B 


B - B 
s m 


B 


and 


m 


n 


rT 


The effect of variation of B and n was studied and all the 

r r 

other parameters were kept constant- The above said parameters 
were varied one by one keeping others constant to see the effect 
of the varied parameters- 

3.11 COMPirr ATI ONS 


The governing equations were solved by explicit method and 
further for solving ordinary differential equation Runge 
K.utta's fourth order method was applied- The whole channel reach 
was divided into 20 cells along length with each ceil of length 
500 m. There was only one cell along transverse direction in 
main channel- For the value of equal to 4, the number of rows 
of cell over flood plains was 4 whereas for = 8 it was S- The 


width of 

each 

cell was 

50 m- the 

value of time 

step was taken 

equal to 

2 sec 

There 

were four 

sets of data as 

shown below. 

Set 

1 : 

B = 4 

n 

o 

1 

11 




r 

r 



Set 

2 : 

B = 4 

n 

= 4.0 




r 

r 



Set 

3: 

B = S 

n 

II 

1 

o 




r 

r 



Set 

4: 

B = 8 

n 

= 4.0 




r 

r 
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In all the above cases unless otherwise specified the values 
of other parameters were taken as ’fol lows: 


dl 

= 50 m 


z 

wop 

18.0 

^f 

= 5 m 


"l = 

0.03 

d 

X 

= 500 m 


9 

9.81 m""/sec 

S 

o 

= O . 0005 


t 

peak 

= 1800 sec 

.^0 

= 5.5 m 


*b = 

3600 sec. 

z 

wo 

= 15.5 m 


^last 

= 28800 sec . 

3.18 EESLfLTS AND DISCUSSIONS 



As 

mentioned earlier. 

four cases 

have been 

studied for 

different combinations of 

B and n - 

The results of the 


r r 


numerical analysis are shown in figures. Figs. 3.5 through 3.8 

show the effect of n^ and on the stage hydrographs at 

different locations at 2 km, 4 km, 6 km and 8 km from u/s 

boundary. Inflow stage hydrograph at u/s boundary has been 

included in each case. Fig. 3.5(a) through Fig. 3.8(a) show the 

stage hydrographs along river cells whereas Fig. 3.5(b) through 

Fig. 3.8(b) show for flood plain cells. Fig. 3.5(c) through Fig. 

3.8(c) show the stage hydrographs at same locations for 

one— dimensional dynamic model developed by liahapatra, K. (1990). 

The plots have been made with normalized time on x— axis and 

normalized depth on y-axis. The figures 3.5 and Fig. 3.7 reveal 

that in the case of two dimensional model there is no subsidence 

in peak depth for the value of n =1.0 where as there is 

r 

considerable subsidence in peak depth in the case of 
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one-dimensional dynamic model. Two— dimensional model has behaved 

like a kinematic model and inflow hydrograph has just shifted to 

next position for the value of n = 1.0. At the same time there 

r 

is not much difference in the shape of stage hydrographs for 

flood plain cells and river cells. Figs. 3.6 and 3.8 reveal that 

there is more subsidence in one-dimensional dynamic model 

compared to two dimensional model even for the value of n = 4.0. 

r 

However, in the cae of n^ = 4.0, ; there is subsidence in 

two-dimensional model too- 

If one compares the effect of , keeping the constant value 

of it is found that there is little effect of B^ on stage 

hydrographs in two-dimensional model- When the value of n is 

r 

equal to 4-0 the time taken to reach the peak depth ' at any 
location is more in two-dimensional model compared to one 
dimensional model at the same location- This is due to ignoring 
dynamic effect. The one— d imensional dynamic model has higher 
wave velocity compared to two-dimensional model without inertial 
ef f ect- 

As it is well known fact that celerity along river is much 
higher than along flood plains, even then there is no difference 
in the shape of stage hydrographs in both of them. This is 
because in river type link equation along longitudinal direction 
bed slope rather than water surface slope has been used. 

Figs. 3-5(a> through Fig- 3.S(a) and Fig- 3.5<b) through ( 
3-S<b) show wavy fluctuations at the end of each stage hydrograph 
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which diminish with time and tend to steady state condition. 
This is due to down stream end ef,fect, because down stream end is 
very near where as in one dimensional model down stream end is 
far away. 

Fig.3.9<a) and 3.9(b> are presented to show the flood peak 
travel time with distance along the river cells and flood plain 

cells, respectively. The plots have been made with normalized 

location on x— axis and normalized peak time on y-axis. The 

(X. 

figures reveal tht n has more effect than B on flood peak 

travel time. It has been also noticed that with larger value of 

the effect is still more. 

Figs. 3.10 show the variation of flood peak depth with 

normalized location. It is clear from Fig. 3.10(a) and Fig. 

3.10(b) for the value of n = 4.0 there is considerable 

r 

subsidence at the downstream locations. In cae of n^ = 1.0 they 
do not behave well and show fluctuation in peak depths in 
two-dimensional kinematic model. However in case of 

one— d imensional dynamic model for all the cases it behaves well. 
In one dimensional dynamic model as shown in Fig. 3.10(c) has 

more effect than n^ on subsidence of peak along channel length. 
When the manning's n for flood plains has been made four times 
the 'n' for main river the peak depth ratio decreases faster i-e- 
the flood plains do play an effective role on flood subsidence 
for both the models- It is also observed that its effect is 
still more with larger value of B . 
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CHAPTER lY 

CONCLUSIONS 


4. 1 CONCLUSIONS 

The various conclusions made are as follows: 

1. The result s of the two-dimensional explicit model for flood 
routine using only continuity principle differ to a great 
extent from the results of the one-dimensional explicit model 
using both continuity equation and equation of motion. Some 
of the results of the two-dimensional model exhibit anomalies 
also. 

2. The two-dimensional model with n =1.0 does not exhibit any 

r 

flood peak subsidence as the flood wave propagates down the 

river. Whereas the model with n = 4.0 exhibits flood peak 

r 

subsidence. The stage hydrographs exhibit oscillations in 
the recession side of the hydrograph. 

3. The time taken to reach peak depth at different locations 

with n = 4.0 in the cae of two-dimensional model is more 
r 

than the time taken in one-dimensional model. 

4. In the two-dimensional model there is not much of a 
difference in the shape of the stage hydrographs of the river 
cells and the flood plain cells at any section for a given 


set of data. 
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5. The variation o-f peak depth with distance along the direction 

of flow exhibits oscillations for n =1.0 where as the 

r 

variation does not exhibit oscillations for n = 4.0. 

r 

6. The variation of time to reach peak depth with distance along 

direction of flow shows that the effect of n is more 

r 

compared to B^. 

4.. a RECOMMENDATIONS FOR FUTURE WORK 

.The anomalies occuring in the present study are due to the 
explicit method used to solve the governing equations. So the 
same case should be studied using implicit method. 



I >. 



rigurc 2.11 Continuity equation of a cell. 
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Fig. 2.3 Arrangement of cells groups 
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Fig. 3. 1 Channel Crosa-eection with flood plains 
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LINK TYPE 


River link 


Weir link 


Side cells 


Intermediate 
Flood plain 
cells 


Next to 
river cells 


River cells 


Next to 
river cells 


Intermediate 
flood plain 
cells 


Side cells 



Fig. 3-3 Topological breakdown of conputatlonal cells. 















































NOR.T I ME(-t/-tpeaK) 

Fig. 3*5 (a) Stage hydrographs In main river for B » A*0, n * 1*0 (2-0) 
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NOR.T I MEKt/t-peaK) 

Fig* 3*5 <b) Stage hydrograjAis In flood plain for = 4*0, = 1*0 (2-D) 



45 



NOR.TIMECt/tpeaK) • model 

Pig. 3* 5(c) Stage hydrograpiis for one — dlraensional/for B b 4*0, n * 1.0 







(oA/A)MJLcd3a'fcdOM 


NOR.T I MECt/tpeaK) 

Fig. 3.6(b) Stage hydrographs flood plain for B s 4*0, n^. * 4*0 ( 2 -^) 
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COA/A)M_Led3a‘fcdOr\l 


NOR.T I MECt/tipeaK) 

6(e) Stage hydrographs for one dlEaensl<Mnal model for B ■ 4*0, n « 4*0 
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(0A/A)M_Lcd3a*td0M 


NOR.T I MEiC-tZ-LpeeiK) 

Fig, 3*7(a) Stage hydrographs in main river for B s» 8.0, n * 1.0 (2-D) 




NOR.T I MECt/tpeaK) 

Fig, 3* 7(b) Stage hydrographs in flood plain for B » 8*0, n ■ 1 »o(2-d) 
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■fcdOM 


NOR.T I MEKt/tpeaK) 

S(a) Stage hydrographs in laain river for ■ 3«0, n^. » 4*0 (2 -d) 
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(oA/A)JHJ-cd3a-fcdOM 


NOR.T I MEICtZ-tpeak) 

Pig. 3.3<c) Stage hydrographs for one dimensional model for . 8.0 
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>H>y3cd-tdOM 


NOR.D I STAMCECx/yo) 

^ig» 3«10(c) Normalized peak depth vs normalized distance £or one 
dimensional model* 
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APPENDIX A 

inr rpnr^OAH FOR FLOOD ROUTINE THROUGH A RIVER WITH FLOOD PLAINS. 
FvrtlCtl METHOD HAS BEEN USED TO SOLVE THE DIFFERENCE EQUATION. 

n»ui^,r-<'»nTA or fourth order has been used to solve 0.0. E. 

THIS FPnCPAM CAN BE RUN DIRECTLY FOR THE AREA HAVING AT LEAST 
tMRF.f flood plain cells EITER SIDE OF RIVER CELL, 
rnp LE«S NUMBER OF FLOOD PLAIN CELLS SLIGHT MODIFICATION IS 
NEEDED WHILE DEFINIG THE IDENTIFIERS 1 . • , IDENT ( J ) . 

NOTAlToNS^USib IN THE PROGRAM ‘ — 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


A 

coo 

CDU 

CCR.CCF 

ox 

OT 

£*i 

or* 1 

fPILM 

rnirr 

c 

IN 

in . irp, 
JN 

TCAM 

PM 

pr 

AL 

OP 

or 

OP 

unfi 

SOf 

I 

TBA'^r 
tPF AP 
H A<*t 
TMA» 

YO 

TOr^ 

iF 

V 

TN 

7U-’‘ 

:uop 

;n 

zu 

?uri 

2UHAV 


CROSS SECTIONAL AREA BETWEEN TWO CELLS. 

COEFFICIENT OF ORIFICE. 

COEFFICIENT OF WEIR. 

UBIGHTINC FACTORS 

LENGTH OF EACH CELL 

TIME INTERVAL FOR COMPUTATIONS. 

WIDTH OF EACH CELl 

TIME INTERVAl FOR WRITING THE RESULTS. 

MAHWINC'S ROUGHNESS OF RIVER CELLS 
MANNING'S ROUGHNESS OF FLOOD PLAIN CELLS. 

ACCElEPAriON DUF TO GRAVITY. 

TOTAL NO OF CELIS IN LONGITUDINAL DIRECTION. 

THE NODES ALONG LENGTH UHERE STAGES ARE BEING CALCULATED. 
TOTAL NO OF CELLS IN TRANSVERSE DIRECTION. 

ROW NO. OF RIVER CELL. 

WETTED PERIMETER FOR RIVER CELLS. 

UEITEO PERIMETER FOR FLOOD PLAIN CELl.S 
Ft.OU ASSOCIATED WITH LEFT FACE OF CELL. 

FLOW ASSOCIATED WITH RIGHT FACE OF CELL, 
now A^iSOCIATEO WITH TOP FACE OF CELL. 

FLOW ASSOCIATED WITH BOTTEM FACE OF CFLL . 

iongitudinai ped si ope for river. I 
longitudinal bed Gl''‘PE FOR FLOOD PLAIN. 

TIME OF COMPUTATION 

BASE OF THE STAGE HYDROGRAFH. 

TIME UPTO PEAK OF THE HYDROGPAPH. 

T1T1E UPTO COMPUTATION IS ProUIRFD 
TIME TO PEACH PEAK STAGE IN THE CELL 
INITIAL DEPTH IN RIVER CELLS. 

INITIAL DEPTH IN FLOOD PLAIN 
HEIGHT OF RIVEP RANK 
DEPTH OVER CELLS AT KNOWN TIME STEP 
depth over cells at NEXT TIME STEP 
STEADY STATE STAGE OF HYDROGRAPH, 

PEAK STAGE OF HYDROGRAPH. 

ETEVATION OF WATER AT KNOWN TIME STEP. 

ELEVATION OF WATER AT NEXT TIME STEP. 

PEAK STAGE IN CELL 


DIMENSION 7B( 30,30>,ZUf 3U , .su i , t i , 

dimension YN( 30 , 30) , IOF.NT (30) , AKH - 
DIMENSION Al 3 ( 30 , 30 > , AK 4 ( 30 , 30 ) , Y T 
DIMENSION y.K TO, 30) # ZWT (30,30) ,ZUe( 
DIMENSION ZUMAXT 10,30 ) , THAX(30,30) 


OPEN (lfMIT'5,rilF» 'RK41 DAT') 

QPPN <t»NI T-EE ,r ILE** 'HURTY OUT M 


p (T A Jll c ♦ ) 

PCAP* 5 . ♦ > 
READ'S. *) 


TN , TN . OX , tn . SO, DT ,PFL T , JCAN 
yO, YF.FRICH.FPICF ,C.rR,CrF 
TPE AK , T BASE , tl AST , 7WO , 2U0P 
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c 


8 


8 

C 


c 

c 

c 


19 


1 $ 

c 

c 

c 

€ 

C 

5 


. * > IC1 , 1C?, IC3, IC4|,ic5,ic6 
Pr . ♦ ) f iPF.NTf ! ) , J- 1 , JN) 


rO»HAt«MM -'.SK.l 5,/, -JN -'^SX^ie,/, *DX -S5X.F6 a,/, 

'PI •* .SV,F6.e,/. -$0 -%5X.F8.6,/, 

•DT -.‘.SX.FS 8,/,'OFLT - * . -<X , F6 . g , / , ' JCAN *.',4X,I8,/, 

•YO .SX,FS.a,/,'YF -'»5X,F6.e,/, 'FRICM -MX,F6. 4,/. 

'FRICF -MXiFe. 4./, *CCR,CCF , 1 X, CF6 . fi, ✓ , ' TPEAK , 1 X , F7 . 0 , / , 

•TFACf •* , 1X,F7. 0,/, 'TtAST - ' , 1 X , F7 . 0 , / , » 2U0 - ' , 3X , F6 . 3, /, 

•ZWOP -^eX.FS 3,/.MCMC2.IC3,IC4,ICS,IC6 -',ax,6I3) 
WR!TF(*.9) ( IDEMK J> . J*| ,JfN) 

FonHAT< * iDCMK j> -'.ax.^ia) 


PH«PL^2. ©♦YF 
FF-Pl 

90H*f t O^FFICH J^SORT tSO) 
FOF-fl 9/FRirF)*SQRT<50) 
+Y01/FH 
YOF - >0 • vr 
«r-| M OF T/FF 
VO^«OM4RH4^0 667 
VOF-ror *RF4*0. 667 
Y ^ 0 9 
T|NT*0 0 


INIMAI CONDITIONS 


00 to 

7U(1,Jt-?W0 

IF nnF?4!t J} NE II THFN 

Y < ! . n • > 0 • V F 

ZP< ! , n-ZU( 1 , J)-Y( 1 , J) 

list 

Yt f , Jl- YO 

zo( t , n«?u( 1 . ji-Yt 1 , j) 

SNDIt 
CON! |NMF 
00 !*•- I-r, IM 
DO »•' 1 t , IN 

ZWMAve 1,11*0 0 

zwt 1 , n -zui ! ' 1 , jr)-Dx*so 
If < lontl! <31 NE. n THEN 
Y< 1 , 11 * YO-YF 
€t. 

YM . J> ■YO 
CMDir 

ZF«I,.n*7W(I,J)-YtI,J) 

CONTINUE 

UM«'.TfAnY CONDITIONS 

UPST*?rAM POUHDRY BY STAGE HYDROCRAPH 


connwjK 

IF< I I.F IPEAKI THEN 
TERM! ’Tt 7U0F-ZU0)/TPEAK 
ZU0UN*ZU0*T£RM14T 
ELSE 

TTPhl-^t ZViOr-ZUOl/ C TEASE- TPEAK) 
;un»iN*7M0P' TERHlfl T-TPEAK) 
ENDir 

IFfT r-l I BASE) 2W0UN*ZU0 
DO *0 1-1 , IN 
2MN< 1 , J )*7U0UN 



o u u 


\ri ir»rMT( n mf n then 

iUn . ?) •» « 7UOUN-7UO)MYO'*Yr ))/YOF 

t\ ^ f 

\Mf » . n- t » 7UOUN'-7UO)*YO)/YO 

Etioir 

go CONTINUK 


UH«iTEAOV CALCUIATJONS AT INTERIOR NODES 


DO r*=‘. l-P.IN -1 
00 tS I-JCAN,JN 
Vy*YI I # J> 

Yl»Yt I -I . j) 

YR-Y( I*t . J) 

YT-Yf I . J«-t ) 

YO-^YI I . J -1 > 

ZTV^^rui I . j) 

Z 0 T** 7 U< I . JO) 

zuo-zun.j-t) 
zzo'-ysn, j) 

7 PT» 7 B< J. J^-t) 

WintNTI J) 

CO TO f 30 . 35 . 40 , 45 ), K 
30 coHTiN»,*r 

I »IVC€lH 0 H 0 T,YY.Yt,YR,YT,ZZy,ZUT,ZBT, 0 X,DL,PM, 50 M,CCR, 

1 CCD 

AWt n . 3 )^ 0 H 0 T 
COlO 

35 continue 

( Ai I NE vt«iVCDHOT,YY,YL, YR, YT,Z 2 W, 2 UB,ZUT,ZZB,DX,DL,PF, 80 F,FR 1 CF, 

1 ct^ . CC^ ) 

AFI M , J t^DHOT 
COTO *“^0 

C At, t. fMf EPMEOIDHOT , YY, YU , YR , YT , YB , ZZU , ZWT , 2 UB,DX, 0 L,PF , SOF ,FRICF, 
1 rcR.tCF) 

AFt ( I , T T^-OMDT 
C f • I 0 

45 ^ 5 Ii 7 'DDf:cnU<DHDT,VY.YL,YR,YB, 2 ZU, 2 WB,OX,DL,PF, 5 nF.FRICF,CCR, 

I ccr ) 

AIM ( f . n^OHDT 
50 CONllN»*r 

YUl, !»' rn. J )^0 S^OT^AKHI.J) 

?U 1 <J . n',U(l,JJ* 0 .S» 0 T»AK 1 «l.J) 
es CONIINl'E 

DO If? I- 3 . 1 N -1 
on If? I'lCAM.JM 
Y M I . .1 I ”• I I I I > 

1 1 ( IM. 11*11 IN, 4 » . 

VY-YH I . J ) 
fi.-ti n -i . Ji 
YP-l 1 1 ! *1 . Jl 
YT’tin . J<n 
YB* Yl I 1 . J -I ) 

Z 7 U -I'M < I , n 
?«t*TU 1 ( 1 . I' I ) 

ZUP* 7 'J 1 « I . I- I I 
rrP'ZPM . J I 
7 ni "fPl ! . IM 1 
P . t Of N I « I I 

cn 10 IIYO.II?, MD.MSI.K 

’ ” riru ' R I VC B t.L I DHBT , Y Y , Yt . YR . YT . ZZU . Z«T . ZB I . OX . DL , PM , SOH , CCR . 

1 CM » 

1 , J ' -OHDT ' 



13 S 


1 


140 


t 


145 


I 


150 

tes 


S 30 

1 


«35 

1 

e 40 

1 


£45 

1 


£50 

££5 


1*^0 


fAi I >|IV rPtV(DHOT. YY 
ccp, (rr > 


, YL , YR , Y T , Z2W, ZUB , ZUT , ZZB , OX , DL , PF , 80F , FR 1 CF 


PY?( I , i)«OHOT 

onto i«:o , 

CONTINUE 

CAll INIHRMEOIOHOT, YY,YL,YR, YT , YD, ZZNi ZUT , ZUB , OX , DL i f F , 80 F « FRXCF « 

CCR , ccr ) 

fiVti I , JWDMOT 

0019 l^’O 

CONY fNUF 

CAll S 1 DEC€LL<OHOT,YY,YL,YR,YB,Z 2 U,ZUB,OX,OL,PF, 80 F,FR 1 CF,CCR, 
CCF> 

AKZC I . JJwOHOT 


CONI iNur 

YF( I . n-m , J>«O.S«DT«AKP(I,J) 

ZU.^< I , .1 1'^ZUi I , J}4 0 .5*DT*AKe( 1, J) 

CONTINUE 
DO P?5 IN-I 

00 frs J-JCAN,JN 
vei 1 . J )-Y< 1 , J> 

Y£< IN. J1-*YI IN. J) 

YY*YP< ! , J> 

YR^Yei l+l . 

YI^YZn , J*1 ) 

Y8^>£< I , J-1 ) 
zzw-'zueii, j) 

2WT'»ZU»U,l’*n . 

ZTm*ZWE( I , J-1 ) 

ZZRvZPC I . J» 

ZBT-ZBI 1 , JM > 

K-TOENT( n 

CO TO «aTO, £35, 240,245) ,K 

COMT iNur 

< ALt R I VCHl L 1 OHOT , YY , Yt , YR , YT , 22U , ZWT , ZBT , OX , DL , PM , 80M , CCR , 

ccr ) 

AK3I f . .n*OHDT 
COTO 

cntniNur 

C A» I YTRIVIDHOT, YY,Yl.,YR, Yt ,ZZW,2UB,2WT^2ZB,DX,DL,PF, BOFiFRlCF, 

CCR. ccr I 

AV II I . JmOHDT 

COTO P*50 

COMT INUE 

CAll INTEPMEDI0M0T,YY, YL, YR, YT, YB,ZZU,ZUT,ZUB,DX,DL,PF,80F,FRICF, 
rCR.Ctr) 

AIMf I , .n->rOHOT 

COTO 

COM! TNUF 

CAll 5IDECEtL(DHDT,YY, YL , YR. YB , Z2W, ZUB, OX , DL , PF , SOF, FRICF, CCR, 

ccr > 

AXK I. JT^OHOT 
GOTO f‘*=0 
COUT INUE 

Y J< I , J »»Y( I , J)+0T*AK3( I, J) 

ZU3I I. T)^ZU(I,J) + DT’^AK3(I,J) 

CONTINUE 
00 3.75 

00 3P5 J^ICAM,JN 
r3( I . J)*Y< t ,.T) 

Y3( IN. JT'^YTIN.J) 

YY- » 31 I , J) 
tL*-Y3( I J» 

YR«y3( 1*^1 ,J) 
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» f I . JM > 

. J ! ) 

I 1 , n 

71IT'*7WH ! , J + 1 ) 

7un*zuM ! . j- n 
7?P-7n» I . J) 

ZBl -7RI ! ,.!♦! ) 

K*ir»rMT<J) 

CO 70 7310, 335, 340. 34S),K , 

330 COWTMItif: 

CAtt RIVrELLCDHDT.YY , YL , YR , YT , 22U, ZUT, 2BT , DX , DU , PM , SOM, CCR , 

I ccr I 

AV.i( T , J )*dHDT 
CO TO 3‘“0 
335 CONTIN’ie 

CAT I NEXTRIVt DHDT, YY, YL, YR,YT, 22U,2UB,2UT,Z7B,DX, DL,PF,SOF,FRICF, 
! CCR,CCP) 

AF4f I . I )«OHOr 
GOTO 3*50 
340 COH7INVF 

TAIL INTFRHFDIDHDI , YY , YL . YR , YT , YB, Z2U, ZUT , ZUB , DX , DL , PF , 80F , FR 1 CF , 

I CrR.CC*^) 

A«4t I . 3 »*0HDT 
GOTO 3*^0 
345 COM7lN'.»E 

CAT. I SID£CEIL<DHDT, YY ,YL,YR. YB,2ZW,2«P,DX,DL,Pr, 50F,FRICF,CCR. 

I ccr ) 

AP4f I , 3>*0H0T 
COlO 

350 COtMINL'E 

385 COMTINUE 

DO 400 I*e. IM-! 

00 4P0 J-JCAN,JW 

07U»f‘A 0*(AKHI,J)+E 0^AK?(I,J)^C.0tAK3( I,JT4AK4(r,jn 
ZUMTt. J)’r;UTI,J)+D2W 
C DEPTH NORMALISED 

IF T 1DENT7 J) NE \ ) THEN 

Y4II I, n'*<ZUNt l,J)-ZB« I.jn/YOF 

£t «r 

YMi I , JI*<2UN( !, JI-ZBt I, JT )/Y0 
CHfUr 

400 COTlTINfTE 

C • - - 

C DOHMeiPEAH BOUNORY CONDITION 

C ..... - 

0*7 I>-JCAN,JM 

r»TM7 IN, J>«ZgM(lN-l , J)-$040X 

ir < infNT < U NE 1 ) THEN 

YfH |H. J)ir(zgN(lN, IN,J))/YOF 

El 

tN< in. J>*( 2UN(IH, J)-Ze( IN,3))/Y0 
ENOff 

55 CONI I N».»r. 

C — - 

C PEt.F|TING THE VALUES 

C - 

DO SO 1-1 . IN 
DO SO 3*JCAN,JN 
ZO< I . n -ZUNl I , J) 

IP ( T »'T*NT I J I NE . f ) THEN 
YT I , J VNl 1 , J)*YOF 

f(T, j7-YNn,J)*YO 
EMf»IT 

SO rONTpn.ip 

C NORMAT ISED PARAMETERS 



TN-^ T /TrpAK 


UR I Tin*, ruir Rg«:in. T? 


IF(T EO !INT) THEN 

URITEt 6S)TN, YNdCa, JCAN> ,YNC 1C3, JCAN), 
YNfir^, jrAN),YN« ICS, JCAMI 
rORHATt 1X,SF0.-O 

URITEf 6?)TN, YN( ICS , JCANf 2 1 , YN< 1 C3, JCAN+S ) , 
YN< IC4 . KAN^e) , YN< 1C5» JCAN42) 
FORMAT(1X,5F8.^) 

TINT-TINT^DELT 

ELSE 

ENOir 

DO 500 I-l , tN-1 
DO 500 J»1 , JM 

IFdUMd.J) GE.ZUMAJCdd)) THEN 
2 UHAi(< I , n*2UN< I , J) 

TMAXf I . J)*T 
ELSE 
ENPir 
tONT INUE 
7-^TfDT 


lEd LE TLAST) THEN 

COTO 5 

ELSE 

URITEt^^^) 

URITEt+.^l ,‘PEAK VALUES AND OCCURANCE TIME* 


UR 1 IE(*.*) 
URITEf < , • ) 
WRITE ( 4 . ♦) 
URITE(t, ♦) 
UR n E C ♦ . ♦ ) 
URITE(t,») 
URlTEI't,'*) 
URITEC t.») 


•AT IC*1C? RIVER NEXIRIVER SIDE ' 
ZUMAXClCa, JCAN) ,TMAX( ICE, JCAN) 
ZUMAXdCa, JCAN+t) .TMAXdCS, JCAN+I ) 
ZUMAXdCe, JCAN’^E) , 7 MAX t ICE, TCAN + E) 
ZUHAXdCE, JCAN+3) ,TMAX( ICE, JCAN+J) 
ZUMAXdCE, JN) ,TMAXdCE, JN) 

» » ^ ... .. ....... 


upnE(*,^) 
WRITEC ) 
UniTF(^,4) 
UP!TE(*,t) 
uniff d.i) 
URIfC ( ^ ) 


'AT IC-ICT RIVER NEXTRIVER SIDE ' 
ZWMAXdC3, ICAN) , TMAX( IC3, JCAN) 
2UMAXdC3, JCAN4 I ) ,TMAX( I f 3, JCAN H ) 
2WMAK(1C3,JCAN4 8) ,TMAX( IC3, ICAM^E) 
7UMAXI IC3, JCANfl) , TMAX( ICI, JCANd) 
2UMAX( IC3. JN) , TMAX< fC3, JN) 


WR n E ( ♦ . * ) • 


UnnF.(*,4) 'AT IC^IC4 RIVER NEXTRIVER SIDE ' 
UP I TE ( * . ♦ ) ZUMAX( ICi , JCAN) , IMAX ( IC4, JCAN ) 
uni TF ( ♦ , ♦ ) 2UMAX( IC<I, JCAN4-n , TMAX( IC4, JCANf 1 ) 
UPITEC ♦ , ♦ ) 2WMAX( 1C4, JCAN^E) , TMAXl IC4, JCAN+E) 
V'PHE ( ♦ d ) ZUMAXC IC4. JCAN43) . TMAy( IC4, JCAN+3) 
UR n E ( ♦ . 4 ) ZUMAX ( lC4,JN),TMAXdC4,JN) 

URI TF ( t , t ) * - 

UPITEd.O 'AT IC-IC5 RIVER NEXTRIVER SIDE ' 
WRITE* ♦.♦ ) ZUMAX ( I C5 , JCAN ) , TMAX ( 1C5, JCAN) 

UR ITS < t , 4 ) ZUMAX* ICS, JCANf I ) , TMAXC ICS, JCAN ♦ 1 ) 
URI TC < t , 4 ) ZUMAX* ICS, JCAH-*E) , TMAX* ICS, JCAN ♦£) 
UR! 7E < ♦ ,4 ) ZUMAX* ICS, JCAN* 3 ), TMAX C ICS, JCAN43) 
URI IFC»,t) ZUMAX* ICS, JN), TMAX dC5,JN) 

URITE<4,4)' 

ENDir 

STOP 

END 


SUtPOUlINE STARIS FROM HERE 


PI'^FP CEI I 

SUPPOUr INE R I VCELL* OHDT , YY , YL , YR , Y T , 22W, ZUT , ZBT , DX , DL , PM , SOM , CCR 



67 


1 ( ^ r ) 

ri,j 0*?QRT(P Rf )*DX 

C?-! 

CJ^Py-'PL 

A*-PL»0 f»(>L*YYI 
ft-A 'PM 

01 *A«P**0 66 7*50M 
A-OLtO 5t<YY4YR) 

R-A ‘ PM 

OR-'A^R^^O 667^80H 
IF<77U CT ZUn THEN 
H-Z7U-ZU1 
U-YT 

row- 0 611 ♦© 0T5^H/U 
CPO*0 6M-0 t75 + U/(HfU) 

IF « f ZZU-ZBT) .CT. I S^YT) THEM 
CDO»0 0 
61 TP 

cn*^**© 0 

ENDir 

OT -CCR* 1 C I ♦COU^Hf^ ! . 5 ♦€?« CDO^U^H^^ 0 ,5) 

El TP 

-77U 
U-77W zei 

rru*o AiMo 07s^H/w 

C00»0 611 '0 . 1 75^W/tH*U) 

IFKZrw-ZBT) IT.0.666 + YT1 THEN 
(00* 0 0 
FI *=C 
rpM-O 0 
FNnir 

ni ^ ( TP# « Cl ♦rDU4H4't1 . 5 + CE + CDO^U+HttO . 5 ) 

FNIMF 

PHOT • I '»l OR -? 0#DT >/C3 

Pt* ttrpii 

END 

C - 

C (Et.t WFvf TO THE RIV6P 

?ltnP0U1 I ME NHXTPIVCDHDT ,YY*YL, YR,YT,2ZW,ZUBiZWT,ZZB»DX, 
PL.pr ,«»or ,FPfcr,ccR,ccr) 

CT-? 0'3 0tSORT<e,0^9 6t)^DX 
t?-*! ?»C» 

Arpt *0 E*<Y11YY) 

Rr A-'rr 

OI^A*R4#0.667^SOF 

E*<YY*YR) 

R-'A PP 

OR'AfP*»0.667**SOF 
IF<77y.LT ZWB) THEN 
H»ZUB-7?»J 
U* Tt 

CHU-0 fill ^0.075*H/U 
(OOrO 611-0. I7S»W/(H+U) 

IF (Y Y IT . 0 . 666*1 ZWB-ZZB)) THEN 

CDO^O 0 

EL$E 

CPU*-0 0 

ENOIF 

OP*U RM C1#CDU*H**I .S+Ca^CDO^UMH^O.E) 

rt It 

M^TTP .UP 
7UP t 7P 
r.pUr A f I I f 0 0 
rno-0 MI-0 175*U/(H*U) 



o o o 


6 8 


77P» LT 0 <.66*YY)THEM 

< nu* 0 0 
ENDIf 

OP--CfP*l Cf *CDW^H^*f s^ce^coo+M^M^^o.si 
ENOir 

A^'DX^O 5<»(YVtYT) 

R-A/DX 

IFC7UT It 2ZU) THEN 
M OPE ^OP T n Z2W- 2UT ) /DL ) /FR I CF 
El ?C 

SI OFF- SORU tZUT“Z2U)/DL)/FRICF 
ENOIF 

OT-CCF ♦A*R'**0. 6S7*SLOPE 

DHOT-<Ol-OR+Q8-QT)/C3 

RFTURN 

END 

C 

c INICPHEOIATE FIOOD BANK CELL 

C 

SUBPOIM INE INTERMEDIDHOT ,YY.YL,YR, YT, YD|2ZU,ZVIT,ZUB, 

I OX .01 . PF .SOr .FRICF,CCR,CCF) 

Cl*? 0/3 0*SORT(8.0^9.8n*OX 
CE-1 

t3*py#nL 

A-Ol *0 5^( YL^YY^ 

Rr A^rM 

I>1 Af 7*SOF 

ArOl «0 SM YY + YR) 
prA#rr 

OP*A^P*tO Af>7*SOF 
A-^pv^n SH YY»YB) 

P" A UY 

ir (isVt I I 7UP> THEN 
?l 0» H-^nPIMZUB-ZZUY/DI. )/FRICF 
FI ?F 

«si OPE* SOPTl (7;W-2UB)/0L )/FR!CP 
ENOir 

OS* ( < F ♦ A^P^tO 6A7*8LOPE 
A*r'> >*y4YT ) 

P 1* A / 0 '< 

If ( rru p ! 7U1 ) YMfN 
et <tr F - ' '»• J < <27U-ZUT )/DH/FRICF 
SI '» 

?| nrr- RIIPM I7UT-ZZU)/0L)/SBICF 

swpjr 

O* 't <r I A*R* *0 , 66?* SLOPE 
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